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La presente invention concerne une methode pour generer un rnaillage 
hybride d'une formation heterogene traversee par une ou plusieurs 
discontinuites geometriques, dans le but de realiser par exemple des 
simulations. 

5 La methode s' applique particulierement a la formation d'un rnaillage 

adapte a un gisement ou reservoir souterrain traverse par un ou plusieurs puits, 
ou par des fractures ou failles, dans le but d'y modeliser des deplacements de 
fluides tels que des hydrocarbures. 

Etat de la technique 

10 La generation de rnaillage est un element crucial pour les simulateurs de 

reservoir de nouvelle generation. Le rnaillage permet de decrire la geometrie 
de la structure geologique etudiee au moyen dune representation en elements 
discrets dans lesquels on effectue la simulation suivant un schema numerique 
approprie. Une meilleure comprehension des phenomenes physiques necessite 

15 de simuler en 3D des ecoulements polyphasiques dans des structures 
geologiques de plus en plus complexes, au voisinage de plusieurs types de 
singularites telles que les stratifications, les failles, les biseaux, les chenaux et 
les puits complexes. Toute cette complexite doit etre prise en compte en 
premier lieu par le rnaillage qui doit restituer aussi fidelement que possible les 

20 informations geologiques dans leur caractere heterogene. 

La modelisation de rnaillage a connu de grands progres ces dernieres 
annees dans d'autres disciplines telles que Taeronautique, la combustion dans 
les moteurs, la mecanique des structures, etc. Cependant les techniques de 
rnaillage utilisees dans les autres domaines ne sont pas transposables telles 
25 quelles dans le monde petrolier, car les contraintes de metier ne sont pas les 
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memes. Par exempie, en simulation de reservoir, les schemas numeriques sont 
construits k partir de volumes de controle afin de mieux respecter la 
conservation de la masse dans le cas d 1 equations de transport de nature 
hyperbolique. Le maillage doit etre de type dit "block-centered" c'est-a-dire 

5 que les noeuds doivent etre situes a l'interieur de chaque couche et les lirnites 
de chaque bloc doivent suivre l'interface entre les couches.. Or si on ne tenait 
pas compte de cette contrainte, on placerait naturellement les noeuds le long 
des failles et le long des lirnites de stratification. Ceci aurait pour consequence 
de faire passer ces interfaces a travers le volume de controle utilise. La 

10 saturation, constante dans le volume de controle, ne pourrait alors respecter la 
discontinuity et les resultats seraient imprecis. II est done necessaire de mettre 
au point de nouvelies techniques mieux adaptees aux besoins petroliers. 

Les maillages cartesiens, couramment employes dans les simulateurs 
commerciaux actuels, sont mal adaptes a la resolution de ces problemes 
15 nouveaux poses par Texploitation des gisements petroliers. Les maillages 
cartesiens, etant a base d* elements parallelepipediques, ne permettent pas de 
representer de telles geometries complexes. 

On connait une methode pour generer des maillages hexaedriques 
structures en 3D de type dit CPG (pour Corner-Point-Geometry) qui respecte 
20 la geometrie des corps. On la trouve decrite dans le brevet FR 2 747 490 (US 5 
844 564) du demandeur et egalement dans la publication suivante : 

- Bennis Ch. et al. « One More Step in Gocad Stratigraphic Grid 
Generation » : Taking into Account Faults and Pinchouts ; SPE 35526, 
Stavanger, 1996. 
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Ce type de maillage offre plus de souplesse que le maillage cartesien car 
il est constitue d'elements hexaedriques quelconques pouvant etre degeneres. 
II respecte rigoureusement les horizons, les failles et permet de representer 
certaines inconformites telles que des biseaux, car sa construction s'appuie sur 

5 ces elements. Toutefois ce type de maillage ne permet pas de resoudre toutes 
les complexites geometriques comme par exemple des maillages radiaux 
circulaires autour de puits complexes. II est possible de realiser separernent le 
maillage du gisement et les maillages autour des puits mais il est difficile de 
representer plusieurs objets dans un meme maillage de gisement de type CPG a 

10 cause des problemes de raccordement lies au caractere structure du maillage. 

Une autre approche est egalement connue par laquelle on genere 
automatiquement des maillages 3D uniquement a base d'elements 
tetraedriques de Delaunay, avec un raffinement radial circulaire autour des 
puits. I/avantage d'une telle approche est qu'elle est entierement automatique 
15 et ne necessite quasiment pas d' intervention de Tutilisateur. Toutefois cette 
methode presente des inconvenients qui rendent les resultats obtenus 
difficilement exploitables : 

- le nombre de mailles est en moyenne cinq fois plus grand que celui d'un 
maillage de type CPG pour une meme structure, ce qui est tres penalisant 

20 pour les calculs de la simulation ; 

- contrairement aux maillages structures qu'on arrive aisement h visualiser, a 
explorer de l'interieur et a modifier localement de fa?on interactive, il est 
tres difficile et parfois impossible de bien controler les maillages 
tetraedriques a cause de leur taille et surtout du fait de leur caractere non 

25 structure. Ceci pose probleme autant pour valider le maillage d'un point de 
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vue geometrique que pour comprendre et valider le resultat d'une 
simulation sur ce type de maillage. 

D'autres approches sont egalement connues permettant de generer des 
maillages, notamment des maillages bases sur des volumes de controle generes 

5 a partir d'une triangulation auxquelles sont associees des techniques 
d'agregation des triangles (ou tetraedres) en quadrangles permettant de reduire 
le nombre de mailles. Bien que des resultats prometteurs aient ete obtenus avec 
ces nouveaux maillages, la representation precise de la complexite geologique 
des reservoirs et des puits reste un sujet de recherche et de developpement. 

10 Malgre leur aspect hybride, ils restent entierement destructures et seraient done 
tres difficiles a gerer et a manipuler en vrai 3D. De plus la prise en compte des 
failles reellement 3D et les puits devies accentuerait grandement cette 
difficulte. 

Par la demande de brevet parallele FR 99/15120 du demandeur, on 
15 connait une methode pour realiser en 3D un maillage hybride sur un milieu 
heterogene qui comporte V utilisation d'un premier maillage structure, de type 
CPG par exemple, et de maillages radiaux structures autour de trajectoires de 
puits ou conduits dessines ou importees dans le modele pour mieux tenir 
compte des contraintes particulieres lies aux ecoulements a proximite de ces 
20 puits. On les combine en inserant les maillages radiaux autour du ou de chaque 
puits dans le maillage global du reservoir apres avoir forme dans celui-ci des 
cavites assez grandes pour y former des maillages de transition non structures. 

Differentes techniques sont connues pour former des maillages non 
structures. On peut les realiser a base de polyedres canoniques (tetraedres, 
25 pentaedres, pyramides, etc.) selon le schema numerique utilise, et Ton connait 
des solutions entierement en 3D applicables a ces types de maillages. 
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La methode selon 1' invention vise a generer en 2.5D, des maillages de 
transition a base de polyedres de type Voronoi* ou analogue pour y appliquer 
des schemas numeriques de type a volume de controle, que Ton ramene a un 
probleme en 2D a condition que toutes les couches des maillages de puits et du 
5 maillage du gisement puissent se projeter verticalement dans un plan 
horizontaLpour former des mailles identiques. 

Un tel maillage doit respecter certaines contraintes : 

- le segment reliant les centres ou sites de deux cellules voisines doit etre 
orthogonal a la face commune aux deux cellules ; 

10 - deux cellules voisines doivent partager entierement la face qui les 
connecte, ce qui revient a dire que par une meme face, une cellule ne peut 
etre connectee qu'a une seule autre cellule ; et 

- les cellules doivent etre convexes. 

Ces contraintes, en 2D, imposent que : 

15 - les aretes des polygones frontieres ne doivent pas etre modifiees (chacune 
doit correspondre a une arete d'une cellule creee, sans subdivision de 
T arete), et que les cellules creees doivent etre convexes ; 

- chaque cellule doit avoir un centre tel que la droite joignant les centres de 
deux cellules adjacentes soit perpendiculaire a l'arete partagee par les deux 

20 cellules ; et 

- les centres des cellules doivent etre situes a Tinterieur de leur cellule. 
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D'autres approches entierement non structurees sont egalement connues 
permettant de generer des maillages, notamment des maillages bases sur des 
volumes de controle generes a partir d'une triangulation auxquelles sont 
associees des techniques d'agregation des triangles (ou tetraedres) en 
5 quadrangles reguliers permettant de reduire le nombre de mailles. 

Les diagrammes de Voronoi sont formes, comme on le sait, en 
determinant autour d'un ensemble de sites P repartis dans un espace, des 
regions formes d' ensembles de points qui sont plus proches de chaque site que 
de tout autre site P de Tensemble. Cette technique de construction est decrite 
10 de fa?on detaillee notamment dans les documents suivants : 

- Aurenhammer, F., 1991 : Voronoi Diagrams : a Survey of Fundamental 
Geometric Data Structure. ACM Comput. Surv., 23, 345-405 

- Fortune, S.,1992 : Voronoi Diagrams and Delaunay Triangulations. Pages 
225-265 of D.Z. Du & F.K. Hwang (eds), Computing in Euclidean 

15 Geometry, 2 nd edn. Lecture Notes Series on Computing, vol. 4. Singapore : 

World Scientific. 

Les cellules d'un diagramme de Voronoi* respectent les contraintes 
suivantes : l'orthogonalite entre les centres et F arete adjacente est respectee et 
les cellules sont convexes. II est done tout naturel de proposer de remplir la 

20 cavite par des cellules de Voronoi, tout en cherchant a respecter les contraintes 
geometriques des bords. Cependant, si Ton cree un diagramme de Voronoi a 
partir des points existants (extremites des aretes des polygones), (cf. Fig. 6) on 
se rend compte que ce diagramme viole les contraintes geometriques, que les 
aretes des polygones ne font pas partie du diagramme et que les aretes sont 

25 coupees 
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Afin d'introduire les aretes des polygones frontieres dans le diagramme, 
il est possible de dupliquer les points existants. De cette fagon, le diagramme 
de Voronoi* est oblige de passer par les aretes. Toutefois, ceci n'est pas encore 
suffisant car si les aretes sont bien dans le diagramme, elles sont encore 
5 divisees (cf. Fig.7). 

Puisqiie les aretes du diagramme de Voronoi correspondent aux 
mediatrices des segments definis par deux sites, on peut essayer de faire passer 
ces mediatrices par les extremites des aretes definies par nos polygones. Pour 
cela, pour un point extremite donne, on choisit les sites du diagramme de 

10 Voronoi* le long des aretes adjacentes a ce point, a egale distance de celui-ci. 
Puisque les points sont dedoubles, on est assure d'obtenir une arete du 
diagramme de Vorono'f sur l'arete de notre polygone. Cependant, si le resultat 
obtenu est correct dans le cas d'un polygone regulier (Fig.8a), ce n'est pas le 
cas pour un polygone non regulier (Fig. 8b). Au moment ou Ton place les 

15 derniers sites, ils ne sont pas forcement a la meme distance du point que les 
premiers que Ton avait places. 

Dans le cas general, l'idee d'utiliser des cellules de Voronoi pour 
remplir la cavite ne semble done pas adequate. 

La methode sdon Tinvention 

20 La methode selon Tinvention permet de simplifier la formation d'un 

maillage hybride adapte a un milieu heterogene traverse par au moins une 
discontinuity geometrique de geometrie connue, dans le but de former un 
modele representatif d'ecoulements de fluides dans ce milieu en accord avec 
un schema numerique defini, la structure du milieu etant connue a priori a 

25 partir de donnees disponibles acquises par des mesures in situ, par des analyses 
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et/ou des interpretations d' images du milieu (des images sismiques par 
exemple, dans le cas d'un gisement). 

La methode comporte la formation d'au moins un premier maillage 
structure (irregulier de type CPG par exemple) pour mailler au moins une 

5 partie du milieu, la formation d'au moins un deuxieme maillage structure pour 
mailler une autre partie du milieu, la formation d'au moins une cavite entre les 
maillages structures englobant chaque discontinuity geometrique et de taille 
suffisante pour que Ton puisse former au moins un maillage non structure 
assurant la transition entre les maillages structures, et la formation de chaque 

10 maillage non structure de transition en utilisant des diagrammes de puissance 
et en imposant aux maillages de transition d'etre conformes au maillage des 
parois de chaque cavite. 

La methode s' applique par exemple a un milieu heterogene ou au moins 
une discontinuite geometrique est un conduit ou puits de geometrie connue 
15 traversant le milieu, auquel cas on forme un premier maillage structure dans le 
milieu en respectant ses discontinuites, on forme un deuxieme maillage 
structure de type radial autour de chaque puits, chaque cavite etant delimitee 
autour de chaque deuxieme maillage radial structure en desactivant des mailles 
de chaque premier maillage structure. 

20 La methode s' applique aussi par exemple a un milieu heterogene ou au 

moins une discontinuite geometrique est une fracture ou une faille traversant le 
milieu, auquel cas on forme un premier maillage structure et un deuxieme 
maillage structure dans des parties du milieu heterogene, de part et d' autre de 
chaque fracture, en respectant leurs discontinuites, chaque cavite formee pour 

25 y inclure un maillage se transition non structure etant formee en desactivant 
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des mailles des premier et deuxieme maillages structures, de part et d' autre de 
chaque faille ou fracture. 

L' utilisation qui est faite ici de la technique des diagrammes de 
puissance est particulierement avantageuse et commode en ce qu'elle permet 
5 de relier entre eux de maniere conforme les maillages structures non reguliers. 

L' invention concerne aussi une methode pour simuler en accord avec un 
schema numerique defini, revolution d'un processus tel que des ecoulements 
de fluides, dans un milieu heterogene traverse par au moins une discontinuite 
geometrique de geometrie connue, en accord avec un schema numerique 

10 defini, la structure du milieu etant connue a priori a partir de donnees 
disponibles acquises par des mesures in situ, par des analyses et/ou des 
interpretations d' images du milieu. Cette simulation est obtenue par la 
formation d'un maillage hybride constitue d'au moins un premier maillage 
structure (irregulier de type CPG par exemple) pour mailler au moins une 

15 partie du milieu, la formation d'au moins un deuxieme maillage structure pour 
mailler une autre partie du milieu, la formation d'au moins une cavite entre les 
maillages structures englobant chaque discontinuite geometrique et de taille 
suffisante pour que Ton puisse former au moins un maillage non structure 
assurant la transition entre les maillages structures, en utilisant des 

20 diagrammes de puissance et en imposant aux maillages de transition d'etre 
conformes au maillage des parois des dites cavites, et en resolvant le schema 
numerique dans le maillage realise sur le milieu pour modeliser le processus. 

Le maillage hybride global est ainsi obtenu par combinaison de 
plusieurs types de maillage : un maillage structure de gisement, un maillage 
25 radial autour de chaque puits, egalement structure, et des maillages non 
structures de transition qui assurent la liaison entre les deux types de maillage 
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precedents. Chacun de ces maillages possede ses propres methodes de 
formation, 

En utilisant un simulateur de gisement d'un type connu, tel que par 
exemple les simulateurs ATHOS™ ou SCORE™, a un gisement muni d'un 
maillage hybride obtenu par la methode, on peut realiser des simulations de 
production. 

Presentation des figures 

D'autres caracteristiques et avantages de la methode selon l'invention, 
apparaitront a la lecture de la description ci-apres d'exemples non limitatifs de 
mise en ceuvre, en se referant aux dessins annexes ou : 

- la Fig.l montre un exemple de maillage structure d'un reservoir faille ; 

- la Fig.2 montre un exemple de maillage radial autour d'un puits vertical ; 

- la Fig.3 montre un exemple de maillage d'un puits horizontal ; 

- la Fig.4 montre en 2.5D, un exemple de reservoir maille ou Ton a menage 
des cavites pour des puits mailles, avant l'etape de creation de maillages 
interstitiels non structures destines a les raccorder les uns aux autres ; 

- la Fig. 5 montre un exemple de maillage en 2D avant creation du maillage 
de transition ; 

- la Fig. 6 montre un exemple de diagramme de Voronoi* qui viole des 
contraintes geometriques ; 

- la Fig.7 montre un exemple de diagramme de Voronoi dont les aretes sont 
divisees ; 



11 



les Fig. 8A, 8B montrent des exempies de diagrammes de Voronoi, 
respectivement dans un cas regulier et dans un cas quelconque avec 
conflit ; 

les Fig.9 a 14 illustrent les definitions et proprietes des diagrammes de 
puissance, rappelees ci-apres ; 

la Fig. 15 montre un exemple de diagramme de puissance obtenu entre deux 
polygones ; 

la Fig. 16 montre un exemple de maillage interstitiel en 2D forme par la 
methode selon F invention, entre un maillage structure de reservoir et des 
maillages structures formes autour d'un puits vertical et horizontal 
traversant le reservoir ; 

- la fig. 17 montre un exemple de maillage global obtenu par la methode ; et 

- la Fig. 18 montre cinq puits pourvus chacun d'un maillage radial, integres 
dans un reservoir maille, par le biais de maillages de transition non 
structures k base de mailles polyedriques quelconques. 

Description detaillee 

La generation du maillage hybride de gisement est effectuee par etapes 
avec addition/sous traction de maillages elementaires de type differents. 

1) Pour representer le gisement global, on utilise un maillage structure 
en i, j, k d'un type connu des specialistes dit CPG, comme decrit dans le brevet 
precite FR 2 747 490 du demandeur. 
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Le reservoir peut etre faille avec glissement d'un bloc par rapport a 
l'autre. Les horizons et les failles majeurs sont d'abord modelises par des 
surfaces continues a partir de donnees issues d'une interpretation damages 
sismiques du sous-sol ou de donnees obtenues durant des forages (marqueurs 

5 aux puits). La structure geologique est ensuite decoupee en blocs failles 
s'appuyant sur ces surfaces. Ces blocs sont mailles individuellement, puis 
rassembles. Pour mailler un bloc, on commence par mailler les surfaces de 
bord puis on peuple l'interieur par interpolation transfinie des maillages 
surfaciques des bords. On applique alors des procedures de relaxation sur les 

10 surfaces de bord et a l'interieur pour harmoniser et regulariser ie maillage. Le 
maillage ainsi obtenu respecte rigoureusement les horizons, les failles et 
permet de representer certaines inconformites teiles que des biseaux. II satisfait 
a toutes les contraintes d'ordre geologique. 

Un tel maillage est effectue en accord par exemple avec la methode 
15 decrite dans le brevet FR 2 747 490 precite. 

2) On dessine de fa<jon synthetique ou on importe les trajectoires de 
puits. On genere alors un maillage radial structure autour de chaque puits pour 
tenir compte des contraintes particulieres liees aux ecoulements a proximite de 
ces puits. 

20 Dans l'exemple represente a la Fig.2, le maillage structure autour d*un 

puits vertical est de type radial circulaire. C'est un maillage de type CPG 
egalement. Sa generation consiste dans un premier temps a echantillonner un 
disque en r, 0 dans le plan horizontal. Le maillage 2D ainsi obtenu est ensuite 
projete verticalement sur les differentes couches du maillage du reservoir. Ici, 
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les i, j, k de la structure matricielle correspondent aux echantillonnages en r, 6 
et z respectivement (cf. Fig.l 1). 

Le maillage autour d'un puits horizontal (Fig. 3) est structure en i, j, k, 
du meme type que celui du gisement a la difference qu'un puits ne peut pas 
5 etre faille. II est egalement obtenu en projetant verticalement sur les differentes 
couches du maillage de reservoir, un maillage 2D appartenant a un plan 
horizontal. 

3) Puis, on insere ce maillage radial autour du ou de chaque puits dans 
ie maillage global du reservoir. Pour cela, on cree dans un premier temps une - 
io cavite dans le maillage du reservoir en y desactivant toutes les mailles en 
contact avec des mailles du puits (Fig.4, 5). L'espace degage entre le maillage 
du gisement et le maillage du puits doit etre suffisant pour realiser 
commodement un maillage de transition. II peut representer par exemple 
environ T equivalent de deux couches de mailles. 

15 4) On genere ensuite un maillage de transition non structure dans cette 

cavite (Fig. 16, 17) pour relier le maillage radial structure autour du puits a 
celui du reservoir qui respecte au mieux les contraintes liees au schema 
numerique. L'utilisateur peut a tout moment revenir desactiver le maillage 
d'un puits en reactivant les mailles de la cavite correspondante dans le 

20 maillage du reservoir. 

Pour realiser le maillage de transition on va utiliser une technique 
connue en soi que Ton trouve decrite par exemple par : 

- Aurenhammer, F.,1987 : Power diagrams: properties, algorithms and 
applications. SIAM J. Compute 16(1), 78-96, 
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dont on va rappeler ci-apres quelques notions pour la clarte de r expose. 
Point pondere 

Soit p un point de St 3 auquel on associe un poids w. On appelle point 
pondere et on note p le couple (p, w) ainsi forme. Un point pondere peut aussi 

etre interprets comme un cercle C de centre p et de rayon sfw . 
Notion de puissance 

On designe par Puis {pt) = Puis {Sx) - \px\' - w la puissance d'un point x 
par rapport au point pondere p ou par rapport au cercle C. JpjcJ" represente ici 

la distance euclidienne entre p et x. Geometriquement, si jc est a Texterieur du 
cercle et si t est le point ^intersection d'une tangente a C passant par x, 
Puis^ x) =\xt\ 2 (cf. Fig. 9). 

Puissance d'un point par rapport a un cercle 
Soit C le cercle de centre c et de rayon r, d'equation : 
C(x) = 0 ou C(x) =xc z - r 
C(x) est la puissance du point x par rapport a C. 
Diagramme de puissance 

Soit S un ensemble de points ponderes. A chaque point p de S, on 
associe une region /?(p), potentiellement vide, definie par : 

K(p)= {re SR 3 1 Puis {p x) < Pttis {Pi ^ Pi e S - {p}} 
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Cet ensemble de regions, ou encore cellules est appele diagramme de 
puissance de S. 

Axe radical de deux cercles 

L'axe radical de deux cercles est la droite formee par les points ayant la 
5 meme puissance vis a vis des deux cercles. En particulier si les deux cercles 
ont une intersection non nulle l'axe radical est la droite passant par 
l'intersection (Fig. 10). 

Utilisation des diagrammes de puissance 

Les diagrammes de puissance se pretent bien a la resolution du 
10 probleme pose car ils ont l'avantage de posseder les proprietes suivantes : les 
cellules qui le composent sont convexes, chaque cellule possede un centre 
(site) et la droite reliant les centres de deux cellules adjacentes est 
perpendiculaire a 1'arete commune aux deux cellules, proprietes qui satisfont 
done a une partie des conditions posees. 

15 On souhaite en outre que les cellules du diagramme de puissance 

respectent les aretes des polygones frontieres. II faut done que les aretes des 
polygones soient des aretes du diagramme de puissance. Pour realiser un tel 
diagramme, on impose comme contrainte que les aretes des polygones soit les 
aretes d'une triangulation de Delaunay couramment utilisee dans le domaine 

20 des maillages. Cette contrainte est facilement realisable car on determine les 
polygones frontieres lors de la creation de la cavite. Une fois cette condition 
etablie, il est possible comme on va le voir, de creer un diagramme de 
puissance respectant les contraintes. 



16 



Diagramme de puissance contraint : respect des aretes donnees 

Soit (s) un ensemble de n aretes dans le plan appartenant a une 
triangulation de Delaunay (T d ) (cf. Fig.l 1). 

Soit (s*) T ensemble des aretes duales de (5) dans le diagramme de 
5 Voronoi" associe a (r rf )(cf. Fig. 12). 

Soit a i et b h i = l..n , les extremites de r arete S i . On considere son arete 
duale S* dans le diagramme de Voronoi associe a (T rf ). Elle definit un 
ensemble de cercles (C ( ) de centres c. pris sur S* et de rayon flc.-aj = 1^^! . 
Tout cercle C ( . de (C.) porte done a ( et 6 ( et, est un cercle vide (Fig. 13). On 
10 considere C] et C; deux cercles de (C, ) . D'apres ce qui precede, a. et b t sont 
sur C) et Cr . a i et b ( sont done les points d'intersection de C) et C; et se 
situent sur leur axe radical. En termes de puissance, on peut done ecrire : 

c 1 l (a,)=cr( ai )=c 1 'fe)=cr(6 1 )=o 

Par construction, quelque soit C } de (c y ), 7=1. .nj*i, Tinterieur de 
15 Cj est vide et ne contient done ni a, , ni b i . On peut encore ecrire : 

vy,c>,)>o 
vy,c>,)>o 

De la, a i appartient a une arete du diagramme de puissance, portee par 
Taxe radical de C) et de C, 2 , b i appartient a une arete du diagramme de 
puissance, portee par l'axe radical de C) et de C; (Fig. 14). 
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Comme par definition, les cellules dun diagramme de puissance sont 
convexes, on en deduit que le segment J= S, fait partie du diagramme de 

puissance. 

Ceci est vrai pour tout S. de (s). On peut done construire un 
5 diagramme de puissance contraint contenant l'ensemble des segments de (s). 
Le degre de liberte pour construire un tel diagramme est grand car la seule 
contrainte est de choisir deux centres par arete duale Sj pour tout S i de (s). 

On remarque que, si S* est reduite a un point, 1'ensemble (C ( ) est reduit 
a un cercle de centre c i = S * . II faut done utiliser dans un tel cas un centre c IV 
io virtuel, infiniment proche du premier pour definir un second cercle. 

Construction du maillage de transition a Taide des diagrammes de 
puissance 

Pour creer un diagramme de puissance, maillage de transition, on 
commence par construire une triangulation de Delaunay des extremites 
15 sommets des aretes des polygones frontieres. On verifie que toutes les aretes 
initiales sont bien des aretes de la triangulation obtenue. Ceci fait, on construit 
le diagramme de VoronoY associe aux sommets des polygones. Ce diagramme 
de Voronoi definit pour chaque arete A des polygones, une arete duale A 1 

Par ia suite, pour chaque arete A, on positionne deux sites du 
20 diagramme de puissance final sur A Puisque les sites doivent etre des points 
ponderes, on associe au site x situe sur A' un poids equivalent a la distance 
||ax|| 2 ou a est une extremite de A. Le diagramme de puissance obtenu a partir 
de ces sites respecte alors les aretes des polygones. La Fig. 16 est un exemple 
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de diagramme de puissance obtenu entre deux polygones avec un seul site par 
arete. 

Une fois ce diagramme obtenu, toutes les contraintes de depart sont 
maintenant respectees sauf a bien verifier que chaque site est a l'interieur de sa 

5 cellule. Par construction, on a positionne les sites sur les aretes duales sans en 
preciser la maniere. On a jusqu'a present rnontre qu'il existait un diagramme de 
puissance correct des lors que les sites etaient choisis sur A \ mais les 
possibilites dans le choix de ces sites sont multiples. II est possible toutefois de 
garantir que si toute arete A intersecte son arete duale A\ il existe au moins un 

10 diagramme de puissance respectant les polygones frontieres et dont les sites 
sont a l'interieur de ieur cellule. 

Resultats 

La Fig. 16 rnontre deux exemples de maillage 2D. Le premier exemple 
correspond a une situation de puits vertical quant a la seconde elle represente 

15 le cas d'un puits horizontal. On peut verifier sur ces exemples que les cellules 
generees pour joindre les maillages des puits et du gisement respectent les 
contraintes enoncees dans le probleme : les cellules sont convexes, les sites 
sont a l'interieur de leur cellule, les aretes des polygones (bords des maillages) 
sont respectees, et la droite joignant les sites de deux cellules adjacentes est 

20 perpendiculaire a i'arete commune aux deux cellules. 

L'etape 2D perrnet dans un premier temps de generer les maillages de 
transition et ainsi de connecter les maillages de puits avec le maillage de 
gisement. Pour obtenir des maillages en volume, on transforme les maillages 
initiaux en maillages 2D, et on applique aux resultats en 2D un traitement 
25 inverse, permettant d' obtenir un maillage global final en 2,5D. 



19 



On a considere le cas ou les deuxiemes maillages structures etaient de 
type radial et formes autour de puits traversant le milieu, avec delimitation de 
cavites autour de chaque deuxieme maillage pour y inclure un maillage de 
transition. II est clair cependant que la methode s'applique au maillage d'un 

5 milieu presentant d'autres types de discontinuity geometrique, tel que par 
exemple un gisement souterrain traverse par des fissures. Dans un tel cas, on 
forme des maillages structures dans des parties du milieu de part et d' autre des 
fissures et Ton delimite des cavites en desactivant des mailles des deux 
maillages structures autour de chaque fissure, assez grandes pour y former des 

10 maillages de transition non structures en utilisant des diagrammes de puissance 
et en imposant aux maillages de transition d'etre conformes au maillage des 
parois de chaque cavite. 
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REVINDICATIONS 



• 1) Methode pour generer un maillage hybride adapte a un milieu 
heterogene traverse par au moins une discontinuite geometrique de geometrie 
connue, dans le but de former un modele representatif d'ecoulements de 
fluides dans ce milieu en accord avec un schema numerique defini, la structure 
du milieu etant connue a priori a partir de donnees disponibles acquises par des 
mesures in situ, par des analyses et/ou des interpretations d' images du milieu, 
comportant la formation d'au moins un premier maillage structure pour mailler 
au moins une partie du milieu, la formation d'au moins un deuxieme maillage 
structure pour mailler une autre partie du milieu, la formation d'au moins une 
cavite entre les maillages structures englobant chaque discontinuite 
geometrique et de taille suffisante pour que Ton puisse former au moins un 
maillage non structure assurant la transition entre les maillages structures, 
caracterisee en ce que Ton forme chaque maillage non structure de transition 
en utilisant des diagrammes de puissance et en imposant aux maillages de 
transition d'etre conformes au maillage des parois de chaque cavite. 

• 2) Methode selon la revendication 1, appliquee a un milieu heterogene 
ou au moins une discontinuite geometrique est un conduit ou puits de 
geometrie connue traversant le milieu, caracterisee en ce que Ton forme un 
premier maillage structure dans le milieu en respectant ses discontinuity, on 
forme un deuxieme maillage structure de type radial autour de chaque puits, 
chaque cavite etant delimitee autour de chaque deuxieme maillage radial 
structure en desactivant des mailles de chaque premier maillage structure. 
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« 3) Methode selon la revendication 1 ou 2, appliquee a un milieu 
heterogene ou au moins une discontinuite geometrique est une fracture ou une 
faille traversant le milieu heterogene, caracterisee en ce que Ton forme un 
premier maillage structure et un deuxieme maillage structure dans des parties 
5 du milieu heterogene, de part et d' autre de chaque fracture, en respectant leurs 
discontinues, chaque cavite formee pour y inclure un maillage se transition 
non structure etant formee en desactivant des mailles des premier et deuxieme 
maillages structures, de part et d' autre de chaque fracture. 

c 4) Methode pour simuler en accord avec un schema numerique defini, 
io 1'evolution d'un processus tel que des ecoulements de fluides, dans un milieu 
heterogene en accord avec un schema numerique defini, la structure du milieu 
etant connue a priori a partir de donnees disponibies acquises par des mesures 
in situ, par des analyses et/ou des interpretations d' images du milieu, 
comportant la formation d'un maillage hybride cons time d'au moins un 
15 premier maillage structure pour mailler au moins une partie du milieu, la 
formation d'au moins un deuxieme maillage structure pour mailler une autre 
partie du milieu, la formation d'au moins une cavite entre les maillages 
structures englobant chaque discontinuite geometrique et de taille suffisante 
pour que Ton puisse former au moins un maillage non structure assurant la 
20 transition entre les maillages structures, caracterisee en ce que : 

- Ton forme les maillages non structures de transition en utilisant des 
diagrammes de puissance et en imposant aux maillages de transition d'etre 
conformes au maillage des parois des dites cavites ; et 

on resout le schema numerique dans le maillage realise sur le milieu pour 
25 modeliser le processus. 
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* 5) Methode selon Tune des revendications precedentes, caracterisee en 
ce que chaque premier maillage structure est un maillage non regulier, de type 
CPG par exemple. 
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